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Abstract. - The generation of entanglement produced by a local potential interaction in a 
bipartite system is investigated. The degree of entanglement is contrasted with the underlying 
classical dynamics for a Rydberg molecule (a charged particle colliding on a kicked top). En- 
tanglement is seen to depend on the structure of classical phase-space rather than on the global 
dynamical regime. As a consequence regular classical dynamics can in certain circumstances be 
associated with higher entanglement generation than chaotic dynamics. In addition quantum 
effects also come into play: for example partial revivals, which are expected to persist in the 
semiclassical limit, affect the long time behaviour of the reduced linear entropy. These results 
suggest that entanglement may not be a pertinent universal signature of chaos. 



Entanglement, i.e. the nonseparability intrinsic to composite systems, is one of the most 
peculiar features of the quantum world. In its most popular form, found for example in the 
original EPR proposal [1], the entanglement is of geometrical nature. However any usual 
potential interaction between two particles can generically lead to entanglement, provided 
several quantum states are accessible to both particles. The study of such dynamical entan- 
glement is particularly interesting for systems which possess a classical counterpart. Indeed as 
is well-known [2], such systems can be investigated with semiclassical tools, allowing to inter- 
pret the quantum dynamics in terms of classical properties. Recent investigations have been 
focusing on the relation between the generation of entanglement and the underlying classical 
dynamics. Initial work in spin-boson systems [3] and in coupled kicked tops [4] suggested 
that chaotic dynamics generate more and faster entanglement. This claim was subsequently 
revised [5,6] and efforts are now being made to derive universal relations ruling the generation 
of entanglement irrespective of system specifics. Different approaches based on perturbation 
expansions [5,7,8], Random Matrix theory [9,10] and semiclassical methods [11,12] have been 
proposed. In particular it was shown and verified [11, 13] that entanglement generated from 
initial Gaussian states averaged over configuration space depends on the global classical dy- 
namical regime. However contrarily to the well established quantum-classical correspondence 
for spectral statistics and large scale fluctuations [14], it is still not clear to what extent an 
analog universal correspondence exists for dynamical entanglement production; for example 
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Fig. 1 - Poincare surfaces of section in the molecular frame, (a) defines the angles 6 and ip for a 
given position of L on the sphere. The core axis is along OZ and N along OX. (e) shows the k=0, 
no coupling situation, valid for all values of T e and To; one sees only the free apparent rotation of 
L around the core angular momentum N, describing circles. On a given circle, the modulus N is 
constant, but it varies from circle to circle, (b)-(d) and (f)-(h): generic case (top row) and resonant 
(T e = To) case (bottom row) for increasing coupling, (b) and (f): k — 0.25, (c) and (g): k — 0.5, (d) 
and (h): k — 10. Note that in the generic case, circles (torii in whole phase space) which correspond 
to low rational values of the returning times of electron and core open resonant islands: values 2 and 4 
appear clearly in (b). Chaos appears along the corresponding unstable manifolds and invades rapidly 
the whole phase space until the last islands of regularity disappear [(d)]. In the resonant case a very 
large couple of resonance islands appears around the ±OZ axis (the top one is visible in (f)-(h)), with 
a separatrix starting from the ±OY axis. 

the dependence of entanglement on non classical features like the type of initial states [7] or 
their spectral width [9] has been put forward. 

In this work we study the generation of entanglement for a real system which possesses an 
unambiguous classical counterpart. This will be done by employing a bipartite system where 
entanglement is produced by a local interaction, viz. the particles interact via a short-range 
scattering potential. Our model enables more stringent tests on the transition to chaos than 
usual models in which this transition depends essentially on a single coupling parameter and 
occurs relatively uniformly on the whole phase space. Indeed we use a system in which a 
second essential parameter, the resonance parameter, produces global changes of the phase 
space picture, changing dramatically the coupling necessary to induce chaos. We will see that 
the underlying classical dynamics is reflected in the generation of entanglement. However this 
dependence involves particular features in the structure of phase-space rather than the global 
dynamical (chaotic or regular) regimes. Quantum effects also come into play, particularly at 
long times. 

The system investigated here is a Rydberg molecule. Simply stated, a Rydberg molecule 
is composed of two parts: a highly excited electron on the one hand, and a compact ionic 
molecular core, containing the nuclei and the tightly bound other electrons on the other. Most 
of the time, the outer electron and the core are spatially well-separated: the core rotates freely 
with angular momentum N and energy oc TV 2 , whereas the outer electron senses a pure 
Coulomb field. Its orbital angular momentum L is fixed in space. Seen in the molecular 
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frame L turns around the N axis (see fig. ^a)). However, the outer electron periodically 
scatters on the molecular core. Classically the electron is kicked by the rotating core. This 
kick results in a change in the direction of L by Aip, the deflection angle of the plane of the 
classical Kepler orbit. TV adjusts accordingly, since the total angular momentum J = L + N 
is conserved. The conservation of L and the cylindrical symmetry of the core implies that 
Ay = k cos 9 where k gives the strength of the kick. A negligibly small kick yields regular 
dynamics. This is visualised in a Poincare surface of section obtained by plotting the position 
of L after each kick (fig. [15]. As k increases, two cases must be distinguished. In the 
generic case [Figs, ^b)-(d)] chaos appears and invades rapidly the whole phase space until 
the last islands of regularity disappear. Things are different in the resonant case. Resonances 
appear when the orbital period of the electron T e and the rotational period of the core T 
are commensurate, so resonances are nongeneric. Chaos is inhibited and finally appears along 
the separatrix [Figs, ^f)-(h)] for much higher values of the coupling k. Such resonances are 
known to modify the spectrum of small Rydberg molecules, giving rise to the stroboscopic 
effect experimentally observed on Na2 [16]. 

In quantum-mechanical terms the collision couples the electron and core dynamics and 
induces phase-shifts in the wavefunction of the outer electron. In the absence of the core- 
electron interaction, the Hamiltonian Hq is separable and the eigenstates are given by the 
product state 

c/> N (E,r)=f L (E-E+,r)®\N) (1) 

where E is the total energy, \N) the state of the core and /l the radial Coulomb function of 
the electron whose energy is ejy = E — E^ [17]. The total Hamiltonian H — Hq + V includes 
the core-electron interaction potential V. The eigenstates are given from scattering theory by 

ij)(E,r) =J2 z n(E)[I + G (E)K]4> N (E,r), (2) 

N 

where Go is the Green's function, Zn are expansion coefficients and K is the scattering matrix, 
whose elements K^'n are related to the transition probability amplitude between states </>jv 
and 0at' due to the collision. From eq. J5J it is apparent that an eigenstate is given by a 
superposition of core states with rotational number N, each core state being associated with 
an outer electron having an energy ejv- The relation between the quantum and classical 
pictures is simpler [15] when expressed in the molecular frame, where K has only diagonal 
elements tan^A- The phase-shifts S\ only depend on the coupling between L and the molecular 
axis. A — L cos 9 is the projection of L on this axis and the classical deflection angle is related 
to the quantum phase shifts by 

*--» £-4 < 3 > 

The first equality is a very general result stemming from the semiclassical approximation to 
scattering theory [18] and is valid to first order in h, whereas the last term follows from the 
functional relation 5\(k) chosen in this work which turns out to be verified for typical diatomic 
molecules. 

To quantify entanglement we will employ the linear entropy S2, as has been done in 
most previous works [3,5-12]. Indeed for bipartite systems S2 behaves as other measures of 
entanglement [20] while being convenient to compute. The linear entropy associated with the 
reduced density matrix describing the electron is given by 



£2 = 1- Tr e/°e> 



(4) 
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Fig. 2 - Average linear entropy for generic (red/gray boxes) and resonant (black dots) states (T e = To), 
shown for different values of k. The rms is also shown as red dashed (generic) and solid (resonant) 
error bars. 

where p e — Tr cj o and as usual Tr c (Tr e ) refers to averaging over the core (electron) degrees 
of freedom. We first determine the degree of entanglement of stationary states. This is done 
by calculating p(E) = \ijj(E)) (ip(E)\ over an energy range for which the classical dynamics 
does not vary appreciably. The values L = 2 and J = 10 adopted in this work correspond to 
experimentally accessible values for diatomic molecules such as Na2. With these values and 
taking into account symmetry requirements, three rotational states are accessible to the core 
N = 8, 10, 12 (classically, 8 < N < 12). Hence in a typical eigenstate three core-electron 
states are entangled. The mean entropy (82(E)) for stationary states is shown in fig. [2] as a 
function of the coupling constant k. The curves for the generic and resonant cases are very 
close despite phase-space being substantially different for most values of k. On average the 
potential interaction entangles at least as efficiently when the underlying classical dynamics 
is mostly regular rather than mainly chaotic. Note the rms is small for the regular and totally 
chaotic regimes and larger in the mixed phase-space situation and that the linear entropy 
saturates in both the generic and resonant cases only for large k: (82(E)) approaches its 
upper limit 2/3 corresponding to a maximally mixed state. 

We next determine the time-dependent generation of entanglement from an initial product 
state. At t = we take the wavefunction of the electron to be radially localised at the outer 
turning point r tp , several thousand atomic units away from the core, and we take the core to 
be in the rotational state |7Vq) , so that 



The mean energy of the radial wavepacket Fi oc is eo ~ E — EJj o and T e is the period of the 
corresponding Kepler orbit. At t rj T e /2 the wavepacket scatters off the core and the wave- 
function is given by a superposition of radial wavepackets moving outward, each associated 
with a core in a given rotational state \N). At t ss T e each of the radial wavepackets has 
reached the outer turning point and starts going back towards the core. The core-electron 
interaction then takes place again at t fa 3T e /2 resulting in further entanglement change. As 
t increases the wavepackets tend to spread radially resulting in a continuous core-electron 
interaction. These features appear in the time dependence of the linear entropy S2(t), ob- 
tained from eq. (jJJ by tracing out the core degrees of freedom from the total density matrix 



For short times S^it) is shown in fig. [31 for Nq = 8 and for 3 different values of the 
coupling k. In the generic case, for small t and not too large k, S2 is seen to increase linearly 
in steps, reflecting the discrete times core-electron interaction, and then appears to saturate 



i/j(t = 0, r) = Fioc(r « r tp ) ® \N ) . 



(5) 
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Fig. 3 - Short time variation of the linear entropy for different values of the coupling strength; t is 
given in units of the Kepler period T e . (a) Generic phase-space, from top to bottom k = 10 (red), 0.5 
(black) and 0.25 (blue). The k — 0.25 (k — 0.5) curve is multiplied by a factor 3 (1.5). (b) Same as 
(a) for the resonant case T e = To. 

(but see below). In the resonant case, the linear entropy also increases on average linearly but 
displays oscillations. Indeed, a semiclassical approximation for the short-time behaviour of the 
purity Tr e p 2 predicts that S^c) k 2 tc- This result, valid only for the collision times tc = 
(2m + l)T e /2, is obtained in the stationary phase approximation by neglecting correlations 
between electron orbits built on different core rotational states. This linear behaviour agrees 
with the prediction of [12]. The k 2 dependence comes from the semiclassical amplitudes, and 
the linear variation in t from the fact that these amplitudes are conserved. The oscillations 
of S*2 observed in the resonant case arise from the coherent projection of the motion of the 
outer electron. It is in this sense a kinematical effect, since no such behaviour is to be found 
in the wavefunction or in the probability flux. By expanding the core energies to first order, 
the purity is seen to depend on terms of the form 

J2 | e -M(;v-JV')/T ] J2 e-^-^ t Z N (E)Z NI (E')(F L (e'(N'))\F L (e(N)))\ 2 : (6) 

NN' e(N)e'(N') 

where ejy = E — E~^ and Fi J (e 1 r) are radial functions on the energy shell (i.e. a combina- 
tion of regular and irregular Coulomb functions with effective phase shifts depending on the 
scattering matrix K). These functions form an overcomplete basis and are not orthogonal 
in e. In the generic case the nondiagonal terms in the second sum are generally small and 
uncorrelated. However at resonance both exponentials are in phase (since e — e' « 2rmr/To) 
and the correlations between nondiagonal terms are enhanced, because the resonance imposes 
a greater similarity of certain corresponding eigenstates (0 [15]. 

The long time behaviour of the entanglement production is shown in fig. ^ In the generic 
case, S2CO is on average very low for k = 0.25 (corresponding to classically regular phase 
space), higher for k — 0.5 (mixed phase-space) and nearly maximal for k = 10 (chaotic phase 
space). Periodic recoherences are visible for small k. Although revivals are expected when 
phase space is regular and the dimensionality of the Hilbert space is low, we point out that 
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Fig. 4 - Same as fig. 3 for longer times; the curves for k — 0.25 and k — 0.5 are not multiplied by any 
factor. The oscillations visible for these curves depend on the revival time of the partial correlation 
functions. For example in (b) for k = 0.5, two scales are visible corresponding to T^ ev (N = 10) and 
the longer Tl ev {N = 12). 



the same observations have been made for considerably larger angular momenta closer to the 
semiclassical limit [19]. The characteristic periods of recoherence involve the revival times Tq™ 
and T£ ev of the core rotation and the electron orbit. The former readily appears by carrying 
the expansion in in eq. JBJ to second order. The latter follows by remarking that the second 
sum in eq. JHJ is the cross-correlation function between parts of the wavefunction belonging to 
different superposition alternatives (i.e. different scattering channels). Alternatively, it can be 
shown [21] that the purity is bounded from below by correlation functions which vary according 
to the revival times. Note that the revivals have no classical counterpart. In the resonant case 
the long-range oscillations can be dramatically enhanced by the commensurate time scales (in 
addition to T = T e , we have Tg e " » 2T£ ev ) and the similarity of the eigenfunctions which 
makes the cross-correlations to be very large. The consequence is that the curves for k values 
corresponding to different underlying classical dynamics may cross. 

These results indicate that the behaviour of the linear entropy is not simply connected to 
the classical regime, but also depends on the structure of phase-space. Indeed the entangle- 
ment production is considerably higher in the resonant case despite classical phase space being 
much more regular (compare Figs, ^c) and (g) with their respective curves in fig. 0J. This 
behaviour can be ascribed to the large values of the cross correlations, which are ultimately 
due to the fact that at resonance L appears to move around the core axis OZ, and therefore 
its projection on TV = OX (and thus the possible values of N 2 ) tend to spread. Classically, 
resonances modify the structure of phase-space by the appearance of a fixed point on the OZ 
core axis, which is visible in the Poincare surfaces of section (fig. [JJ. In quantum mechanical 
terms N is then not well-defined and thus the wavefunctions are mixed in \N) even for small 
kick strengths. In the generic case however these mixings can only take place for large values 
of the coupling, corresponding to chaotic dynamics. 

To conclude we have investigated the generation of entanglement induced by a potential 
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interaction in a bipartite system, a Rydberg molecule. By and large chaotic classical dynamics 
is associated with more and faster entanglement because the scattering particles explore larger 
portions of phase-space leading quantum-mechanically to superpositions. But chaos is not 
necessarily the most efficient way to ensure that the relevant part of phase-space leading to 
superpositions (in our system the rate of inelastic scattering) is explored with a large classical 
amplitude: this is why a dynamical effect such as resonances generate higher entanglements 
despite classical phase-space being mostly regular. On the other hand it has been shown 
[11, 13] that for weak coupling the purity does depend on the global dynamical regime when 
a given type of initial state is averaged over all the possible initial positions. But it is also 
known that the average entanglement is sensitive to the form of the initial state [7]. From a 
general standpoint we would like to point out that the universal relations between quantum 
fluctuations and average classical phase-space properties are known to be spoiled by nongeneric 
families of periodic orbits [22] that are expected to proliferate in systems (such as Rydberg 
systems) that are strongly coupled by potential scattering [23] which is here the generator of 
entanglement. 

* * * 

We are grateful to T. H. Seligman (UNAM, Cuernavaca, Mexico) for fruitful discussions. 
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